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Abstract 

The elastic proton-deuteron backward reaction is analyzed within a covariant 
approach based on the Bethe-Salpeter equation with realistic meson-exchange in- 
teraction. Lorentz boost and other relativistic effects in the cross section and spin 
correlation observables, like tensor analyzing power and polarization transfer etc., 
are investigated in explicit form. Results of numerical calculations for a complete 
set of polarization observables are presented. 
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1 Introduction 



Presently a wideranging program of precisely investigating the structure of the lightest 
nuclei is under consideration. There are new proposals to study the polarization character- 
istics of the deuteron using both hadronic and electromagnetic probes (cf. [0, 3). Be- 
sides the goal of checking fundamental results of quantum chromodynamics (for instance, 
the study of the evolution of the Gerasimov-Drell-Hearn sum rule [^) the envisaged 
experiments focus on a complete reconstruction of the amplitude of the corresponding 
process |I|, ^, |[ and an overall investigation of the nuclear momentum distribution 

The simplest reactions with hadron probes are processes of forward or backward scatter- 
ing of protons off the deuteron. The extensive experimental study of these reactions has 
started a decade ago in Dubna and Saclay (cf. 0, ||, H]) and is planned to be continued 
in the nearest future at COSY |0. One may classify this type of reactions as inclusive 
break-up processes, exclusive quasi-elastic scattering and elastic processes. A common 
feature of these processes is that in a collinear geometry, the measured momenta of the 
fragments are directly connected with the argument of the deuteron wave function in the 
momentum space, supposed the reaction mechanism is dominated by the one-nucleon ex- 
change. In such a way a direct experimental investigation of the momentum distribution 
within the deuteron in a large interval of internal momenta seems to be accessible. By 
using polarized particles one may investigate as well different aspects of spin-orbit inter- 
action in the deuteron and obtain hints on the role of non-nucleon degrees of freedom 
in the deuteron wave function, like A isobars, A^iV excitations and so on. An encourag- 
ing fact here is that the extracted momentum distributions from different reactions with 
electromagnetic and hadron probes are rather similar, and therefore a realization of ex- 
perimental programs at different facilities may provide a quite complete information on 
the internal structure of the deuteron. 

Nowadays the elastic proton-deuteron (p-D) backward scattering with both polarized 
protons and deuterons receives a renewed interest [0]. A distinguished peculiarity of 
this process is that within the impulse approximation the cross section is proportional 
to the fourth power of the deuteron wave function, contrary to the break-up and quasi- 
elastic reactions which are proportional to the second power of the wave function. This 
makes the processes of elastic scattering much more sensitive to the theoretically assumed 
mechanisms, and even a slight modification of the deuteron wave function may result in 
significant deviations from the calculated cross section. However, as pointed out by Vasan 
[Pm], Frankfurt and Strikman |TI|] and Karmanov the polarization observables, like 



tensor analyzing power and polarization transfer, are exactly as those obtained for the 
break-up and quasi-elastic scattering in the non-relativistic limit. Beyond both the non- 
relativistic limit and the impulse approximation the polarization observables differ for 
different processes. Hence, a combined analysis of data on polarization characteristics 
from the above mentioned three processes will constrain the basic reaction mechanism 
and the role of non-nucleon degrees of freedom and relativistic effects in the deuteron. 
Another peculiarity of elastic backward or forward pD reactions is that the amplitude 
of the processes is determined by only four complex helicity amplitudes, and a complete 
reconstruction of these amplitudes seems possible in one experimental set-up. For this it 
is sufficient to measure 10 independent observables as proposed in refs. |l3l. Certainly 
this does not mean that the realization of such an exhausting experiment will determine 
entirely the deuteron structure; only within the non-relativistic impulse approximation 
the cross section is directly related to the deuteron wave function. 

First measurements of polarization observables, such as the tensor analyzing power T20 
and polarization transfer k, have been performed in Dubna |14|, |15| and Saclay ^. 



Theoretically the elastic pD scattering has been studied by many authors |jT6|, [T^, |T9 
|20[| . It has been shown that the cross section can not be satisfactorily described within 
the non-relativistic impulse approximation and that other mechanisms, e.g., described 
by meson-exchange triangle diagrams are important. Besides the importance of 



other mechanisms, the role of relativistic corrections within the impulse approximation has 
been studied by several authors already some time ago (see |18] and further references 
therein) within the Bethe-Salpeter (BS) formalism. The unpolarized cross section and 
the tensor analyzing power have been numerically computed in a fully covariant way and 
a comparison with polarization data, available at this time, has been made. However, in 
view of the present experimental situation and future proposals |^, |^, a detailed 
covariant investigation of the role of relativistic corrections, such as Lorentz boost effects 
and contributions of negative-energy waves etc. is still lacking. In the present paper an 
attempt is presented to fill this gap. 

We focus here on a detailed study of the elastic pD amplitude within the BS approach 
by using the numerical solution obtained with a realistic one-boson exchange interaction 
p2| . It is known that one of the unpleasant features within the BS formalism is the 
cumber someness of the final expressions for the calculated observables and difficulties in 
their physical interpretation (cf. refs. [17, ^ 2^|). In this paper we try to avoid this 



problem and present our results in a form as simple as possible. For this sake we separate 
the contributions of the positive-energy waves and identify them in the non-relativistic 



limit. The contributions of the Lorentz boost effects and the negative-energy waves in 
the deuteron are then calculated in leading order. That means, in the fully covariant 
results we keep the first orders of negative-energy P waves and the leading-order terms of 
a Taylor expansion. In this way we are able to separate and to investigate in explicit form 
the contribution of the Lorentz boost effects and negative-energy waves to the amplitude, 
cross section and polarization observables as well. Results of numerical calculations for a 
complete set of observables are also presented. 

Our paper is organized as follows: In section II the kinematics of the process is described, 
and the covariant amplitude is derived in details in section III. Since the covariant ex- 
pressions within the BS approach are rather lengthy and since in the procedure of the 
arrangement of results in parts containing non-relativistic formulae and boost effects and 
relativistic corrections separately, it is very important to provide a clear definition of all 
the relevant variables. In sections IV - VI a complete set of polarization observables is 
defined in terms of BS wave functions. The non-relativistic limit, Lorentz boost effects 
and relativistic corrections are studied analytically and numerically for the cross section 
and for a class of polarization observables. A comparison with available experimental 
data is also made. In section VII an interpretation of the relativistic corrections in terms 
of non-relativistic meson-exchange like contributions is performed. To do so we solve 
the BS equation for the negative-energy P waves in the one-iteration approximation and 
express explicitly the ^P^' and ^Pi^ waves via the non-relativistic S and D waves of 
the deuteron. The obtained result for the elastic amplitude is found in a form being very 
similar to amplitudes computed in the non-relativistic picture when estimating the role of 
NN pair currents in electromagnetic processes. Some cumbersome expressions and useful 
formulae are collected in the Appendices A and B. 



2 Kinematics 

We consider the elastic backward scattering reaction of the type 

p + D=p'{e = 180°) + D'. (1) 

The differential cross section of the reaction (|1|) in the center of mass system (c.m.s.) of 
colliding particles reads 

where s is the Mandelstam variable denoting the total energy squared in the c.m.s., and 
Ai is the invariant amplitude of the process. In the case of backward scattering the cross 



section eq. @ depends only on one kinematical variable which usually is chosen as s. 
Other variables can be expressed via s by using energy conservation. For instance, the 
Mandelstam variable u is u = (Mj — m'^Y/s, while the c.m.s. momentum is = —t/4 
etc. Here Md and m stand for the deuteron and nucleon masses, respectively. Then in 
the c.m.s we define the relevant kinematical variables as follows: the four-momenta of 
particles read 

D = iE,p), p=ie,-p), D' = {E,~p), p'=(e,p), (3) 

and the polarization four- vectors of the deuteron with polarization indices M and M' can 
be written as 

iM - Cm + PmJeTWY ^ ^ 

where C, C' are the three-polarization vectors in the rest frame of the deuteron, 

C+i = (-l,z,0)/v^, C_i = (l,^,0)/v^, Co = (0,0,1). (6) 

The proton spinors are normalized as u{p)u{p) = 2m with 

/ 



s) = -y/m + e 



_^ <P',s') = V^e\ ^' I, (7) 



where Xs denotes the usual two-dimensional Pauli spinor. 

In what follows we shall widely exploit as variable the momentum of the outgoing proton, 
Piab, in the rest frame of the incoming deuteron, which we here define as laboratory 
system. The relation between c.m.s. and laboratory system is simply expressed by |p| = 
2mPiab/ a/u. 

Now we proceed with an analysis of the general properties of the invariant amplitude 
M.. In principle, the amplitude for the elastic fermion- vector-boson scattering has 
been studied in detail and is well known (see, for instance refs. 0, 13, ^), nevertheless 



for the sake of completeness we present here some of the most important characteristics 
oiM. 

The process of the elastic pD scattering is determined by 12 independent partial ampli- 



tudes |2^. However, in case of forward or backward scattering, due to the conservation 
of the total helicity of colliding particles, only four amplitudes remain independent, and 
these four amplitudes determine all the possible polarization observables of the process. 



There are many ways of representing these four amphtudes. In order to emphasize ex- 
phcitly the transition between initial and final states with fixed helicities it is convenient 
to represent A4 in the c.m.s. in a two-dimensional spin space for the proton spinors and 
three-dimensional space for the deuteron spin characteristics. In this case the manifest 
covariance of the amplitude is lost. However the analysis and final formulae become much 
simpler and transparent. Moreover, by making use of eqs. (§) - (0) all the polarization 
characteristics of the reaction may be expressed via the corresponding quantities evalu- 
ated in the deuteron rest frame. This gives another advantage of such an analysis, namely 
it allows for a straightforward non-relativistic limit, in particular avoiding the problem of 
boosting or not the non-relativistic polarization vectors from the rest frame to the c.m.s. 
It is worth emphasizing that in such a representation of the invariant amplitude, in spite 
of the fact that it is not explicitly covariant, its form is the most general one and valid in 
both the relativistic picture and the non-relativistic limit as well. 

In this paper we keep our notation as close as possible to the one used in refs. |13 . 
Hence the total amplitude is written in the form 

M = xt'^Xs (8) 

with 

^ = ^ i^AdU + 13 {n^M){n^U + ^C{cT■ [^,, x |+ ,]) + zV (o-n)(n ■ [^,, x ^+ ,]), (9) 

where n is a unit vector parallel to the beam direction; A, B, C and V are the partial 
amplitudes of the pD elastic scattering process depending on the initial energy. Then the 
cross section @ is determined by Tr(jF+jF). In the following we suppress the subscripts 
M and M' of the polarization vectors ^, bearing in mind that in computing observables 
the summation over these indices results in the completeness relation for C, or in the 
polarization density matrix of the deuteron, which reads in covariant form |26| 



Ee^.a/^(-.. + ^), (10) 

+ {-\ ((W^aJ,,(W^aJ^ + (W^aJ,,(W^aJ^) (11) 

where {W\)^^ = ie^^^xD^ / M^; Sd is the spin vector, and Qd stands for the alignment 
tensor of the deuteron. /x, i/, A ■ ■ ■ are Lorentz indices, and we use the metric g^^, with 
signature —2. 



In the three-dimensional representation for the deuteron polarization and two-dimensional 
Pauli matrices for the nucleon states the corresponding density matrices can be cast in 
the simple form 

p, = i(/+(^TPp)), (12) 

p''' = l{s^,-l^e-''S^-2Q'^f), (13) 

where Pp is the proton polarization three- vector, Sd and Qd are spin and tensor polar- 
ization operators (actually 3x3 matrices) of the deuteron. 



3 The one-nucleon exchange mechanism 

We investigate here the relativistic one-nucleon exchange defined by the diagram depicted 
in fig. 1. Using the kinematics shown in fig. 1 and working with the Mandelstam tech- 



nique 1 27], the one-nucleon exchange contribution to the elastic amplitude within the BS 



formalism is 

M = u{p') T{D, q) S2f{D\ q') u{p). (14) 

T{D, q) denotes the BS vertex function of the deuteron; 5*2 = 1/ [d/2 — q + rnj is the 
nucleon propagator, and T = 7or"'"7o. We use the abbreviation D = D'^j^ when contacting 
a four-vector with Dirac matrices. The momenta q and q' are fixed by the conditions 
D/2 + q = p' and D'/2 + q' = p. The vertex function r{D,q) is the solution of the BS 
equation for the deuteron bound state. The BS equation and consequently its solution 
T{D, q) are sixteen-component objects in the spinor space. To solve the BS equation and 
to compute observables within the BS formalism one usually represents the vertex function 
r(D, g) as a 4 X 4 matrix and utilizes a decomposition of r{D, q) over a complete set of 
matrices in the sixteen-dimensional spinor space. As mentioned in ref. [Q, the choice of 
the representation of the matrices depends on the special attacked problem. Actually in 
some calculations it is convenient to combine two representations, namely the complete 
set of Dirac matrices (to perform explicit numerical calculations) and the so-called p spin 
classification of the partial vertices |^ (to have a more transparent physical interpretation 



of the obtained results; for details consult ref. p3, p9|). In the present paper we use mainly 



the p spin classification, however the numerical calculations have been performed in terms 



of the solution of the BS equation obtained in the Dirac basis pT] , p^ , |30| . Notice that 



the solution of the BS equation has been obtained in the deuteron rest frame, whereas 



considering the amplitude ([IID it is seen that at least one vertex T{D, q) is defined in a 
system where the deuteron is moving. Obviously in this case one may explicitly boost 
the vertex r{D,q) from the rest frame to the c.m.s. and analyze the amplitude ([T^) in 
terms of a BS solution at rest and Lorentz boost effects separately. Our experience shows 
that in this way one obtains rather cumbersome expressions and a relatively simple and 
physically meaningful analysis of results is straitened. Therefore, we represent the BS 
vertex function in a covariant form with eight invariant scalar functions (see Appendix 
A). Observing that in the process depicted in fig. 1, in each vertex r{D,q) one nucleon 
is on the mass shell, the part contributing to the amplitude (0) may be written in a 
form exactly coinciding with the one used, for instance, by Gross or Keister and Tjon 

ma-' 

r(A,) = |ft.? + feM), + |,,,| + ,,,M|M^l±^, ,15) 

m mm 
where hi are invariant scalar functions depending on the invariants and Dq. Then after 

summation over the spins in eq. (0) one gets the cross section in a fully covariant form in 

terms of the relativistic BS solutions T{D, q) as 
da 1 1 



dVt 647r2s 6 



Tr 



{p' + m)T{D,q)S2T{D',q'){p + m)T{D\q')S2T{D ,q) , (16) 



where, when replacing T[D, q) by eq. (15) and summing over the deuteron spins, one 
should make use of eqs. (|T0|) or (|ll|) in dependence of the initial and final states of the 
particles in a concrete measurement. 

Eq. (jT^) together with eqs. (^) and ([IT]) completely determine all characteristics of the 
process. When computing the trace in eq. ( |16D one obtains a fully covariant relativistic 
expression of all observables. In our calculations we use a suitable algebraic formula 
manipulation code which, within the representation of the solution of the BS equation 
in the form (plsf) , delivers the covariant, but rather cumbersome results (for examples 



cf. ref. |r7|). These results have been tested by evaluating the non- relativistic limits for 
the cross section and polarization observables. However, further investigations relying 
on these lengthy analytical expressions seem to be almost impossible. Therefore, for an 
explicit study of the infiuence of the Lorentz boost effects and other relativistic corrections 
we shall investigate different aspects of the amplitude M. in eq. (0) instead of the cross 
section (pIB]). 

By substituting eq. (^) into eq. (0) and making use of the Gordon identity, 
u(p', s')[d{p — m) + ip' — m)a]u(p, s) = 0, one can represent the amplitude in the form 

6 

M = J2Riu{p')Riu{p) (17) 

1=1 



with six invariant scalar functions Ri and six covariant spin structures defined as 

Ri = it R2 = ——, R3 = i— + i'—, R, = i{D-f)i', (18) 

mm mm 

R^ = Pl_t^^D-p'), R, = i(D-p')^ + iD-pX^, (19) 
mm mm 

R, = -{2h,hl + hhl) - ^ , (20) 

R2 = -{KK + - 775 — ^ — J, 21 

m \^ ~ P ) " ^ 



~ 1 Tnh -\h* 

R3 = -ihK + h2K + hhl) - — i^^, (22) 



6 Ms .9o^ 

[D — p ) — m-^ m^^ 

^5 = ^t: 2 + (24) 



{D — p'Y ~ ""^^ ^ 

-Re = 77:; 2 ~^ 2~- (25) 

At first glance there seems to be a contradiction between eqs. ([T8|) and (|19]) and the 
general expression (P), namely instead of four amplitudes our result contains six different 
structures. However it is straightforward to prove that the six covariant spin structures 
in eqs. (|l^) and (p!9D in the coUinear kinematics reduce to exactly four independent forms. 
For instance, taking into account that, in case of forward or backward elastic scattering, 
the expression {D — p') has no spatial component in the c.m.s. the three structures R2, R5 
and i?6 are equivalent and determine the amplitude B in eq. (^ (see below) . The structure 
may be cast into the form of eq. (^ by exploiting Dirac's matrix algebra, 

ii' = m - ^ui = (eeo + 70(^17) + ^7075(^27), (26) 

Pi = (eof + a"), V^ = eappU'p; a, (3, p= 1,2,3. (27) 

Then it is seen that among the six amplitudes eqs. ([T8|) and ([T9|) only four are indepen- 
dent in collinear kinematics. The correspondence between eq. (P) and eqs. ([TSD and ( p!9| ) 
becomes obvious if the former are written in the c.m.s. In this case. 



u{p',s)Riu{p,s) = x1> 



i^^')2e-zcT{^x^')^2m (28) 



Md{Md + E){m + e) Ma Md 



X. 



Si 



u{p',s') R2u{p, s) 
u{p',s') Riu{p, s) 



,2 



Xs 



E-e 



Mdm 

X X. 



(29) 
(30) 



MdMd 



Xs 



u{p', s') R4^u{p, s) 



X 



u{p',s') R5u{p, s) 
u{p', s') Rqu{p, s) 



{Md + E-m- ef 
Md{Md + E){m + e) 



2{E-t 



+ t(T{i X $! 



-(^02m + 2(^p)(p,^ x|') 
2{e^ -m^ -Ee) 



Md 



Xs) 




(31) 
(32) 
(33) 



MdMd 

It is worth stressing again that, in spite of these six covariant spin structures Ri - ■ ■ Rq 
being written in c.m.s. have lost their exphcit covariance, they still determine the covariant 
amplitude Ai of the process. The corresponding invariant scalar functions Ri defined by 
eqs. (pO[) - (p5|) may be computed in any reference frame. Since the numerical solutions p2| 



of the BS equation have been obtained in the deuteron rest frame, we also express Ri in 
this system. When one nucleon is on the mass shell, i.e. Md/2 + qo = E'^, the invariant 
functions hi are of the form (see Appendix A) 



1 1 



2^'+2M,Pz,,^" P 



{2E^ - Md)g. 



55 



47r/i2 



y2(m + E;) 

\/^mE'^ 

-95, 



m{m + 2Ep 



2MdPi 



lab 



m 



^91 + 



{E'p + 2m)m' 
2MdP?. 



-97, 



(34) 
(35) 

(36) 
(37) 



V2Md{m + EX^ ' 2MdP;ik ' 2MdPiab 

where Qi are the BS vertex functions in the deuteron rest system and all the kinematical 
variables in eqs. (|3^) - ( p7|) should be evaluated in this system. 



4 Observables 



Having determined the amplitude by eqs. ([T7| ) - ( ^71) one may define various polarization 
characteristics of the process. Employing the notation used in refs. |2^ we define 



the set of all possible polarization observables for the non-covariant amplitude @ by 

_ Tr(^aAPH^+avPHO 
rt\,H^\',H' Tr (J'J^+) ' ^ ' 

where the subscripts A and H (A' and H') refer to the polarization characteristics of the 
initial (final) proton and deuteron respectively; ax is the Pauli matrix, and T> stands for 
a set of 3 X 3 operators defining the deuteron polarization. Note that the introduced sub- 
scripts may appear as either single index or double indices in dependence on the reaction 
conditions. For instance, 0,0 0,0 means a process with unpolarized particles, while 
0, NN — i> 0, NN means the tensor-tensor polarization of the initial and final deuterons 
parallel to the normal of the reaction plane direction. 

At this point it is worth mentioning that the numerical solution of the BS equation 
has been obtained in the Euclidean space-time with imaginary time component go of the 
relative momentum q. In the process under consideration go is fixed and real. Hence, one 
needs either a numerical procedure for an analytical continuation of the amplitudes to the 
real relative energy axis (cf. [0) or another recipe |3^ for using the numerical solutions 



in this case. 



We rely on the analysis of the BS partial vertices performed in ref. where the 
dependence of S and D wave vertices upon the relative energy is shown to be smooth, 
contrary to the amplitudes which display a strong dependence on go. Therefore, in our 
calculations, we can replace, at moderate values of go, the S and D vertices by their 
values at go = with good accuracy. The P vertices can be expanded into Taylor series 
around go = up to a desired order in po/m. Then the corresponding derivatives can be 
computed numerically along the imaginary axis since they are analytical functions of go 



Finally, to cast our formulae in a more familiar form, known from non-relativistic cal- 



culations, we introduce the notion of BS wave functions 17 



^s{\Pia,\) = A/" ".r"r;" > ^^(ip...i) = ^If^' (39) 

^p.(|P.a.|)=Ar^I%l^, (40) 



^?i(0,| 


Plab\) 




-Md ' 




P/ab ) 



where A/" = l/ATTy/2MD- Then the cross section (ITB|) and the amplitude (0) may be com- 



puted in terms of positive and negative-energy wave functions, "^3,0 and '^Ps,Pr respec- 
tively. The BS wave functions "^3,0 are intimately related to the famous non-relativistic 
deuteron wave functions u{Piab), w{Piab), and at small values of Piab they practically coin- 



cide p4|. Therefore, their contribution to the amplitudes and cross section is henceforth 



referred to as the non-relativistic result. The parts containing the negative-energy waves 
\l/pg are of a purely relativistic origin and consequently they manifest genuine relativistic 



correction effects. The weights of these wave functions are quite small [0, Q and we shall 
neglect all the terms proportional to "^j^^ p^; only interferences between "^3,0 and ^I/p5,P7 
are kept in our results. Besides the mentioned relativistic effects there is another source 
of relativistic corrections, namely the so-called Lorentz boost effects stemming from the 
transformation of the BS wave functions from the c.m.s. to the deuteron rest frame. 



5 Lorentz boost effects 

A straightforward way to study effects of the Lorentz boost is to compare the non- 
relativistic results with those obtained in the BS formalism by equating to zero all contri- 
butions from the negative-energy waves leaving only the contribution of the waves \E's,d- 
In the non-relativistic case, the cross section and the spin amplitudes take a simple form 

EM 

'^ = 3{u\q)+w\q))\ (41) 

Anr = [uiq) + ^) , (42) 
Bnr = -lw{q) (2V2uiq) - wig)) , (43) 



Cnr = (n(g) + (n(g) - V2w{q)) , (44) 

T^NR = ^w{q) (u{q) + ^) , (45) 

where in eqs. (^) - (^) the relative momentum q and appropriate sign conventions are 
introduced in the definition of the non-relativistic wave functions u{q) and w{q). Note that 
in non-relativistic calculations there are ambiguities in treating the internal momentum 
q. It is not clear in what frame of reference the argument of the non-relativistic wave 
functions u{q) and w{q) should be evaluated: is q to be computed in the three-body 



center of mass frame as proposed in ref. [0] or should one use the deuteron rest frame 



[pO|] in defining the deuteron wave function? At small momenta these alternatives become 
identical, however a difference occurs already at intermediate energies Piab ~ 0.2 ■ ■ ■ 0.3 
GeV/c. In the covariant BS approach this problem is solved by using invariant amplitudes 
hi{q^,Dq) and by taking into account the boost effects. 
Substituting eqs. (^) - (pTl) into eqs. ([T6| ) and (^) and expanding the result into 



Taylor series around P^^^/2m^ and keeping the leading terms we obtain the contribution 
of positive-energy BS waves in the form 

dao _ 12m2 . ^,^2.p ( . , PL , ^QPL . SSP^, \ 



Ao = IGnmPL [<fs{Piat) - ) ma,), (47) 

Bo = lenmP^l'^DiPlab) (2V2^s{Plab) + ^D{Plab)) C{Plab), (48) 



2 



2 



Co = 167rmP2, Us{Plab) - ^^^) {"^siPlab) + V2^D{Plab)) C{Pi 



lab) 



(49) 



Do = - 167rmP,l, ^^^DiPlab) [^^siPlab) - "^""^ J /:(P^a.), (50) 

where the Lorentz boost effects are represented by C{Piab) defined as 

^(«-) = (i + jt?+*?? + -)- 

It is seen that the results within the BS approach recover the non-relativistic formulae in 
the leading order in P;^^/2m^ and receive additional corrections from the Lorentz boost 
effects. At small values of Piab these corrections are negligible but at moderate values 
of Piab they become important and may give up to 30 - 40% contributions in the non- 
relativistic cross section. Note that within the one-nucleon exchange mechanism, Piat is 
kinematically restricted so that the Taylor expansion in eqs. ( ^61) - (|5lD is justified since 
P;^^/2m^ < 1 in the whole range of the initial energy ^/s. 

The elastic cross section for a process with unpolarized particles evaluated within the 
BS formalism is shown in fig. 2. The dashed line is the contribution of only positive- 
energy BS functions and (to be compared with the dotted line which represents 
the computation of the cross section in the non-relativistic limit with Bonn potential |35[|). 



The long-dashed curve represents the corrections from pure Lorentz boost effects. The 
solid line is the result of full BS calculations by exploiting the numerical solutions 



2^ of the BS equation with a realistic kernel with vr, lu, p, a, rj, 6 exchanges. Experimental 
data are taken from refs. |2^, ^ . The relativistic effects coming from the negative-energy 
waves are much smaller and they are not displayed here. From fig. 2 it becomes clear that 
the Lorentz boost effects become essential at Piab > 0.5 GeV/c. This is an understandable 
effects, since it is expected that the boost corrections should increase with increasing initial 
energy. Remind that the region Piab > 0.5 GeV/ c corresponds already to rather high initial 
energies, say T^in ~ 6 GeV. (The kinematics of the process is so that at ^/s oo the 
momentum of the detected slow proton becomes Piab 0.75 m.) 



A comparison with experimental data shows that the one-nucleon exchange mechanism 
alone does not describe satisfactorily the cross section and that other mechanism should 



be considered [|T^, One should keep in mind, however, that measurements of the cross 
section in the strict backward direction are rather difficult and many experimental data 
are presented as extrapolations of da/dQ obtained in nearly backward direction to the 
exact backward angle |[16| , In view of the very strong angular dependence of the cross 
section the data from different groups differ noticeably (see fig. 2 in ref. and fig. 27 



in ref. [0). This, together with the above mentioned sensitivity of the results to the 
chosen deuteron wave function, generates uncertainties in a detailed comparisons with 
data. It is not our aim here to improve the agreement with the data, but we intend to 
proceed with our methodological study within the well defined framework of the impulse 
approximation . 

The Lorentz boost does not affect at all the polarization characteristics defined by 
eq. (|38[) , as it should be. For instance, a direct calculation of the cross section eq. ( |T6|) 
with the density matrix (|ll|) results in the known non-relativistic formulae for the tensor 
analyzing power T20 = —Ho^NN^ofi/ and the polarization transfer k = 37io,L^L,o/2 

rpNR _ 

~ V2 mPlab) + niPlab) ' ^ ' 

V2 niPlab) + niPlab) ■ ^ ' 



Observe that in the non-relativistic limit the polarization characteristics eqs. ( |52D and 
(|53D in the elastic pD elastic scattering exactly coincide with ones in the reactions of 
inclusive [|, || or exclusive deuteron break-up processes. In the relativistic case such 
simple relations among polarization observables do not hold. From this one may conclude 
that a combined analysis of data obtained within different processes would allow for an 
estimate of the role of relativistic corrections in the deuteron wave function. 



6 Relativistic corrections 

In this section we present results of numerical calculations of the contribution to the 
polarization observables coming from the negative-energy P waves. We are going to 
investigate the tensor analyzing power and polarization transfer. 

There are two factors of smallness in computing relativistic corrections: the term 
Pf^f^/2im? and terms proportional to the negative-energy P waves in the BS amplitude. In 
estimating corrections to eq. (^) we expand the result of the trace operation in eq. ( |T6|) 

1 A 



and our results of calculating eqs. (^8]) - (|33| ) into Taylor series around P;^^/2m^ and 
keep only the leading terms relative to negative-energy waves and to P;^^/2m^. The final 
results for the relativistic effects reads for the cross section 



da d(jQ 

dn ^ In ^ ' 



(54) 



Sa 



24\/6m3p, 



2y2P, 



3m^ 



P7 



and for the tensor analyzing power 



20 



20 



20, 



(55) 



ST, 



20 



2^3 (^s + V2^d) {^s - -^d/V^ 



m 



and for the polarization transfer k 



Pi 



^P5 + 



lab 



+ ■ ■ 



.NR 



+ 5k, 



(56) 



It should be stressed that a comparison of the relativistic corrections (|55|) and (|56D with 



those obtained for the deuteron break-up reaction p9| ^ demonstrates that in the elastic 
proton-deuteron scattering the polarization observables coincide with the ones in the 
deuteron quasi-elastic and break-up processes only in the non-relativistic limit. Therefore, 
the data from these processes all together (cf. [0]) can determine the magnitude of the 
relativistic corrections and further constrain the deviation of the mechanism of these 
processes from the simple one-nucleon exchange picture. 

In fig. 3 results of calculations of the tensor analyzing power T20 in the elastic pD back- 
ward reaction are presented. The long-dashed curve is the contribution of only positive- 
energy BS waves eq. (^2]), the dotted line represents the pure relativistic corrections in 
eq. (|55|) , while the solid line is the total result within the BS formalism. Experimental 
data (open and full circles) are from refs. |]^, |n|, |T^. For the sake of completeness we 
present some experimental data for the deuteron break-up processes (triangles) and 
results of computation of T20 within the minimal relativization scheme ||38|, ^ with the 



Paris deuteron wave function In fig. 3 it is seen that the theoretical calculations 

predict a change of sign in the T20 whereas the experimental data from both deuteron 



break-up and elastic scattering processes shows that T20 remains negative in the whole 
interval of measured values of Piab- However, calculations with the positive-energy BS 
wave functions (long-dashed line) together with the relativistic corrections (dotted line) 
result in a better description of data. 

In fig. 4 similar calculations are presented for the polarization transfer n. The dashed 
line is the contribution of the positive-energy BS wave functions eq. (^3|), while relativistic 
corrections eq. ( p6D are presented by the dotted line; the solid line is the total BS result. 
As in the previous results the relativistic corrections in k are small for small and moderate 
values of Piab, but become essential at higher values of the initial energies. 



7 Other spin observables 

As mentioned in ref. 0, one of the goals of the future experiments is a direct reconstruc- 
tion of the four complex amplitudes (^. For this one needs to measure 7 independent 
observables, however as seen from the above formulae, all polarization observables are bi- 
linear combinations of the amplitudes (P) so that the necessary number of measurements 
at given energy increases. A full set of polarisation observables for complete measurement 
has been proposed in refs. |^, |13]. It is found that 10 spin observables, i.e. the two of the 



first order, like cross section and tensor analyzing power TCo^nn^o,o, and 8 spin correla- 
tions of the second order, for instance Ho,nn^o,nn , 'Ho,nn-.o,ss , 'Ho,n^o,ls , Ho^N-^Nfi , 
Hn,n^o,o , 'Ho,LS^N,o , ^o,L5^o,o and 'Ho,ln^o,ln could provide a complete analysis of 
the spin amplitudes (^. 

In the previous section the results of relativistic calculations of the tensor analyzing 
power and polarization transfer from the initial deuteron to the final proton have been 
presented. In this section additional calculations of spin-correlation observables of the 
second order are performed. We consider here the proton vector-vector transfer, the 
deuteron vector-vector and tensor-tensor transfer coefficients. They read explicitly 

, 4 
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where the first hnes in each of eqs. (0) - (^Uj) display the non-relativistic limit, while the 
second lines are the corresponding parts of purely relativistic corrections. 

Figs. 5-7 show these spin-correlation observables according to eqs. (|57|) - (^) calculated 
within the BS formalism; they are depicted as solid lines. The dashed lines show the 
contribution of the positive-energy BS waves (i.e., the non-relativistic limit), while the 
dotted lines are the relativistic corrections due to the contribution of P waves in the 
deuteron. It is seen that the relativistic effects for the proton-proton transfer coefficients 
are negligible (see fig. 5) while for the deuteron-deuteron correlations these effects are 
essential at Piab > 0.5 GeV/c. 

It is interesting to notice that there are observables which in the non-relativistic limit 
are exactly zero and therefore consist of relativistic corrections only. For instance the 
tensor-tensor transfer coefficient 'Ho^ln^o,ln is predicted to vanish in the non-relativistic 
case, while within the BS formalism one gets 

Tr (^^+) no,LN^o,LN = 9 (A^ + ReAB - C^) (60) 

2 



where all the contributions from higher orders in Pf^^/rn^ have been neglected. In spite 
of the small value of this tensor-tensor correlation (see fig. 8) the measurements of such 
observables may directly quantify the importance of relativistic P waves in the deuteron. 

As a conclusion of this section we emphasize that for some of the computed spin ob- 
servables the relativistic effects are not too important, whereas for a specific class of 
observables (like the tensor-tensor correlations) admixtures of P waves in the deuteron 



may result in important corrections, for instance the coefficient (pO| ) differs from zero only 
due to relativistic effects in the deuteron. 
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8 One-iteration approximation 



The relativistic corrections in eqs. (^) - (pOf ) are governed by negative-energy P wave 
states in the deuteron. This can be considered as a hint that admixtures of P waves within 
the BS approach are related to relativistic corrections by taking into account meson- 
exchange currents and NN pair production diagrams |^ in the non-relativistic picture. 



To establish a correspondence between our results and the mentioned non-relativistic 
calculations we estimate the contribution of the relativistic corrections by computing 
the P wave vertices in the so-called "one-iteration approximation" The gist of this 



approximation is as follows: In solving the BS equation by an iteration procedure one puts 
as zeroth iteration the exact solution of the Schrodinger equation for 5* and D vertices and 
zero for other waves; then the P vertices are found by one iteration of the BS equation. 
Our experience in solving numerically the BS equation shows that it converges rapidly 
for relatively small momenta < 1 GeV/c. That means when utilizing the exact non- 
relativistic solutions, after one iteration the resulting P waves are not too far from the 
full solution. 

To obtain analytical expressions for the negative-energy waves we proceed with the 
mixed BS equation in the sense that the BS vertices are expressed via the BS amplitudes 
as 

where G{k) and are the BS vertex and amplitude respectively, A^g^ denotes the 
meson-nucleon coupling constant and jmes the meson-nucleon coupling vertex (for scalar, 
pseudo-scalar or vector couplings). 

Then using the decomposition of G{k) and $(p) in the complete set of spin angular 
matrices Fq (for details consult ref. [|2^) the BS equation for partial vertices Qa and am- 
plitudes $o (here a accounts for the p spin indices of the corresponding partial amplitude) 
may be written in the form 

ga{ko, \k\) = i J Jc^dnkj^—^^—^Tl [r+(fco, -k)lmesT (3{j)0,V)lmes\ $/3(P0, \p\) ■ 

^ ^ (62) 

Using the standard decomposition of the meson propagator over generalized Legendre 
functions Qi{z), where z = (|p p + |A;p + fi'^ — (po — ^o)^)/|p ||^| , one gets 

g^iko, \k\) = -XLs J ^^^W^fs{\k\, \p\) $,(po, IpD, (63) 

Wat3{k,p)= J2 ^li'^ii / c?^fcC?f^p>lm(p)>lm(A;)Tr [r+(A;o, -A;)7mesr/3(po,p)7m, 
1- 27i\p\\k\ J 
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These expressions are still the exact BS equation in the ladder approximation. Further 
we assume: 

(i) in the first approximation the negative-energy waves are exactly zero, i.e., in eqs. ( |63D 
remain only $/3 = S^^ and D^^, 

(ii) in the interaction kernel Wa/s and in the vertex functions g'^^ip) we neglect the 
dependence on the relative energy, i.e., Wa/3{k,p) ^ M/a/3(|/i;|, \p\) and g^^{p) — g^^iO,p), 

(iii) the negative-energy waves are obtained by only one iteration of eq. (§3|). 
For the pseudo-scalar isovector exchange we get 
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where TV = -^/{2\p\\k\Ek), Ek = VWT 



We perform further calculations in the coordinate space: 
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where \I's(^) and \1'd('") are the deuteron wave functions in the coordinate space. 

Then the result for the function "^p^^ with a BS kernel with pseudo-scalar one-boson 
exchange reads 
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(70) 



where u{r) and w{r) are the non-relativistic deuteron wave functions in the coordinate 
representation, and g^ ^ 14.5 is the pion-nucleon coupling constant. The normalization 
factors are = V2 (1) and = -1 {V2) for (^pj waves. 

With this definition of the negative-energy waves one may estimate the origin of the 
relativistic corrections computed within the non-relativistic limit as additional contribu- 
tion to the impulse approximation diagrams, such as meson exchange currents and NN 
pair production currents. As an example we compute within the one-iteration approxi- 
mation the amplitude A which turns out to have the simple form of a negative-energy 
wave contribution 



A = An 



32V67rP3jM/ 
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Substituting ([70|) into the expression for the amphtude eq. (|7lD , the relativistic corrections 
in the one-iteration approximation become 
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MM 

oo 
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jdr^{l + fir)i,irPiab) [v^ n(r) - w{r)] (^^s - j ' (^2) 

which is similar to expressions obtained in non-relativistic evaluations of the so-called 
"catastrophic" and pair production diagrams in electro-disintegration processes of the 
deuteron |^ and which is also similar to results of computation of the triangle diagrams 
usually considered in the elastic pD processes [|l9i |20[| . In our case these corrections 
may be represented as diagrams with meson exchange due to anti-nucleon degrees of 
freedom in the BS equation, as depicted in fig. 9. The remaining amplitudes B, C, V 
and, consequently, the cross section (Q) and all the polarization observables (^) receive 
analogous corrections. From this it becomes clear that generic relativistic calculations, 
even in impulse approximation, contain already to some extent specific meson-exchange 
diagrams, i.e., pair production currents, and one should pay attention on the problem of 
double counting when computing relativistic corrections beyond the spectator mechanism. 



9 Summary 

In summary, we present an explicit analysis of various relativistic effects in elastic back- 
ward scattering of protons off deuterons within the Bethe-Salpeter formalism with a real- 
istic interaction kernel. To have a well defined framework for our methodological investi- 
gations we rely here on the impulse approximation. This allows to identify and investigate 
separately the contributions of the positive-energy waves, Lorentz boost corrections and 
relativistic effects due to negative-energy waves. Particular attention is paid to the com- 
putation of the four spin amplitudes of the process within the Bethe-Salpeter approach. 
By writing these amplitudes in the center of mass system in a non-covariant form a direct 
correspondence between our approach and the general phenomenological analysis of the 
process is found. In such a way a suitable representation of the polarization observables 
and a straightforward investigation of the non-relativistic limit are achieved. 

Numerical estimates of the Lorentz boost and other relativistic effects in the cross section 
and selected polarization observables, at kinematical conditions of ongoing and forthcom- 
ing experiments 0, are presented. It is found that in a complete set of polarization 
observables, proposed for a reconstruction of the amplitude, relativistic corrections either 

on 



may be negligible for a certain class of spin-correlations or play a crucial role for other 
observables. It is shown that the one-nucleon exchange mechanism alone does not give the 
predominant contribution in these reactions and future experiments must clarify effects 
beyond the impulse approximation. 
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Appendix A 

The covariant expression of the vertex function r{D,q) for the BS equation with one 
particle on mass shell is rather known and may be found, for instance in refs. [T7_, 
Here we present the covariant solution T{D, q) for the full BS equation when both particles 
are off the mass shell 

T{D, q) = [hi^ + h2 H [hs^ + /i4 + 

mm m 

^ {qO^ D/2-q + m , I)/2 + q-m ^^ ^ + m 

mm m mm 

The eight invariant scalars hi{Dq,q'^) are connected with the corresponding spin-orbit 
momentum vertex functions gi{po, Piab), which are numerically determined in the deuteron 
rest frame, via 

= ,.1 (2^ - 2po + Mrf)(Mrf + 2po + 2El)gi + 
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where po = {.D<i)/M(i and E^ = y P^^b + The spin-orbit momentum vertex functions 
are defined within the p spin classification as gi = 5*++, g2 = S , gs = g^ = D , 

95 ='P^-, 96 ='P-^, 97 ='P^-, 98 ='P-^- 
Eqs. (|A 2|) - (|A 9|) may be written in a more compact form, however the present expres- 



sions are more informative and easily understood by the reader. 



10 Appendix B 

In obtaining eqs. (^Sf ) - (|33|) the following relations are useful 

(|a X b|)(|c X d|) = (ac)(bd) - (ad)(bc), (B2) 
-!o-(p X ?)(p?') + X ?')(pS) = X - i(crp)(p. C X «'), (B3) 

(<.a)(W) = + X «')^ + f g - 'M.i X Oj^^j^. (B5) 
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Figure 1: The one-nucleon exchange graph for the reaction p+ D = p'{Q = 180°) + D'. 




Figure 2: The spin averaged differential cross section da/dQ for the elastic proton- 
deuteron backward scattering in the c.m.s. as a function of the momentum of the detected 
proton in the laboratory system. Dashed line: contribution of the positive-energy BS 
waves, long-dashed line: contribution of the Lorentz-boost effects eq. (|^), solid line: full 
BS calculations, dotted line: results of calculations within the non-relativistic limit with 



the Bonn potential wave function. Experimental data from |20 




Figure 3: The deuteron tensor analyzing power T20 for the elastic proton-deuteron 
backward scattering. Long-dashed line: contribution of the positive energy BS waves 
eq. (|5^) , dotted line: purely relativistic corrections computed by eq. (^), solid hne: 
results of computation within the BS approach eq. (|55|) , short-dashed line: results of 
computation within the minimal relativization scheme with Paris potential wave 
function. Experimental data: circles - elastic backward scattering ||, |1^, [1^, triangles 
- T20 measured in the deuteron break-up reaction 0. 
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Figure 4: The polarization transfer n for the elastic proton-deuteron backward scattering. 
Notation as in fig. 3. 
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Figure 5: The vector- vector polarization transfer coefficient from the initial proton to 
the final proton. Dashed line: contribution of the positive-energy BS waves (i.e., the 
non-relativistic limit), dotted line: relativistic corrections, solid line: full BS results via 
eq. (13). 
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Figure 6: The vector-vector polarization transfer coefficient from tlie initial deuteron to 
the final deuteron. Notation as in fig. 5 
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Figure 7: The tensor-tensor polarization transfer coefficient from tlie initial deuteron to 
the final deuteron. Notation as in fig. 5. 
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Figure 8: The relativistic corrections for the tensor-tensor polarization transfer coefficient 
from the initial deuteron to the final deuteron 'Hq,ln^q,ln defined by eq. (0). In the non- 
relativistic limit 'Ho,ln^o,ln vanishes. 
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Figure 9: A possible intermediate mechanism already included in the one-nucleon ex- 
change diagram in the BS approach. 
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